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Resumen 


En la mecánica de fluidos existen diversos fenómenos que exhiben un 
comportamiento histerético y con cambios súbitos no sólo de orden cuantitativo, 
sino también cualitativos. Tales fenómenos pueden ser representados y 
explicados mediante la teoría de la catástrofe, que permite explicar este tipo de 
respuesta observada en algunos sistemas físicos. En este artículo se demuestra 
analíticamente que el flujo alrededor de cilindros sumergidos en un flujo 
constituye una catástrofe del tipo en pliegue. Ampliando este resultado, ahora 
con un enfoque semi-empírico, se puede demostrar que el flujo alrededor de 
cilindros y esferas con rugosidad se puede representar mediante una geometría 
de catástrofe tipo cúspide. Ambos resultados respaldan la hipótesis de que el 
flujo alrededor de objetos sumergidos en líquidos se puede modelar mediante 
geometrías de catástrofe y explicar así la ocurrencia de bifurcación, bimodalidad, 
saltos bruscos e histéresis. 


Palabras clave: aplicaciones teoría de la catástrofe, teoría de la catástrofe, flujo 
alrededor de objetos sumergidos, histéresis en mecánica de fluidos, mecánica de 
fluidos, catástrofe cúspide. 


Introducción 


En 1972, René Thom publicó el libro titulado 
Stabilité Structurelle et Morphogénese (Thom, 
1974). En esta obra expuso por primera vez, de 
manera sistemática, los principios de la teoría 
de la catástrofe. Desde entonces, dicha teoría 
ha adquirido importancia como una técnica útil 
para modelar fenómenos de comportamiento 
discontinuo e histerético. La teoría de la 
catástrofe se ha aplicado principalmente en las 
ciencias sociales y biológicas, en menor grado en 
las ciencias exactas y con menor frecuencia en 
la mecánica de fluidos e ingeniería. No obstante, 
con cada vez mayor frecuencia se amplían las 
aplicaciones en diversos campos, por ejemplo 
para la operación de presas (Kuman y Chauhan, 
2009) y en la modelación de flujos de debris (Cui 
et al., 1996). 

La teoría de la catástrofe puede aplicarse a 
la modelación de fenómenos reales de manera 
cualitativa O cuantitativa. En particular, la 


teoría de la catástrofe tiene sus más amplias 
aplicaciones de manera cualitativa. Desde un 
punto de vista conceptual, predominantemente 
descriptivo, valiéndose de las características de 
las superficies de equilibrio de las funciones de 
catástrofe, más que a un método matemático 
riguroso, la teoría de la catástrofe se utiliza de 
manera cualitativa cuando las ecuaciones que 
gobiernan al fenómeno físico en estudio son 
desconocidas o no se dispone de soluciones. 
En particular, debido a la falta de soluciones 
analíticas de las ecuaciones de Navier-Stokes, 
las aplicaciones analíticas de la teoría de la 
catástrofe en mecánica de fluidos son escasas y 
se extiende más un enfoque semi-empírico, en 
el que se emplean de manera descriptiva las 
diversas geometrías de catástrofe. 

El propósito de este documento es 
evaluar, mediante el empleo de la teoría de 
la catástrofe, de manera tanto cuantitativa 
como cualitativa, la hipótesis de que el flujo 
alrededor de objetos sumergidos se puede 
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describir mediante modelos o superficies de 
catástrofe. Para ello será necesario describir 
algunos de los conceptos fundamentales de 
la teoría de la catástrofe. 


La teoría de la catástrofe y la mecánica de 
fluidos 


Conceptos fundamentales de la teoría de la 
catástrofe 


Se dice que la función f.R">R tiene un punto 
crítico en p si: 


of 


OX 
Lp 


of 


OX 
21 


0 (D 


Los puntos críticos de una función son 
topológicamente importantes debido a que 
definen su forma cualitativa y, por tanto, el 
comportamiento general del fenómeno real que 
representan. 

El tipo de puntos críticos de una función 
puede establecerse mediante el determinante de 
la matriz hessiana, conocida también como la 
matriz de estabilidad: 


Y 
Hf= ¡=1,2,3,..., 
f | dx) | j n (2) 
Si el factor determinante de esta matriz no es 
nulo: 
det Hf+0 (3) 
en los puntos críticos p, entonces éstos 


Geométricamente 
hablando, los puntos críticos no degenerados 


no son degenerados. 


son máximos, mínimos o puntos silla. Al 
aplicar el teorema de Morse (1931), este tipo 
de funciones se puede reducir a la forma 
cuadrática: 


f=Y xi (4) 
i=1 


Donde A,son los autovalores (o eigenvalores) 
de la matriz de estabilidad y se conocen como 
funciones de Morse. 

Cuando la matriz de estabilidad es singular 
en un punto crítico u, esto es, cuando: 


det=Hf/,=0 (5) 


el punto crítico se denomina degenerado. 
Bajo estas condiciones, si no se aplica la teoría 
de la catástrofe, se puede decir poco sobre 
el comportamiento de la función alrededor 
del punto crítico. El teorema de Morse no es 
aplicable y, en principio, no existe ninguna 
forma canónica a la cual se pueda reducir la 
función. 

Los puntos críticos degenerados son de 
mayor interés debido a que están relacionados 
con inestabilidades de la función. Se dice que 
una función f.R">R es estructuralmente estable 
si para cada perturbación lo suficientemente 
pequeña p:R">R, la función perturbada f+p 
mantiene el mismo número y tipo de puntos 
críticos que la función original f. Esto es, si existe 
un difeomorfismo que transfiera fa f+p. 

En un sentido geométrico, esto significa que 
la función perturbada f+p mantiene el mismo 
número y tipo de puntos críticos que la función 
original f esto es, que es cualitativamente igual. 

Las funciones con puntos críticos no 
degenerados, o funciones de Morse, son 
estructuralmente estables. Las funciones con 
puntos críticos degenerados son estructural- 
mente inestables. 

Cuando la función representa un fenómeno 
físico, el concepto de inestabilidad estructural 
implica que el fenómeno real modelado por la 
función se puede modificar significativamente 
con una pequeña perturbación. Este tipo de 
comportamiento es conocido desde hace mucho 
tiempo en mecánica de fluidos. Birkhoff (1950) 
citó ciertos casos de flujos con esta propiedad y 
estableció que existe la posibilidad, en mecánica 
de fluidos, de que “causas arbitrariamente 
pequeñas puedan producir efectos finitos”. 
De hecho, en un fenómeno estructuralmente 
inestable, vemos no sólo cambios cuantitativos, 
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sino también cambios cualitativos, bajo 
pequeñas perturbaciones. 

Generalmente, una función tiene en su 
mayor parte puntos críticos no degenerados, 
dependiendo de los valores de sus parámetros y 
sólo en algunos intervalos tiene puntos críticos 
degenerados. 

Para una familia de funciones f.R"xRISR, 
con n variables, llamadas variables de estado, 
y con parámetros k, llamados parámetros de 
control, Thom estableció la existencia de formas 
canónicas si 1<2 y k<5. Estas formas canónicas 
se conocen como funciones de catástrofe. El 
cuadro 1 muestra las funciones de catástrofe 
elementales. 


La condición: 
Vf(x;,0;)=0 (6) 


define la 
estado de cada función de catástrofe. Estas 


superficie de equilibrio o de 


superficies son de especial importancia, ya 
que el comportamiento de histéresis y de 
discontinuidad, entre otras peculiaridades, se 
puede establecer y modelar en ellas. 

La figura 1 muestra la geometría de la 
catástrofe en cúspide, donde las principales 
características de los fenómenos de catástrofe 
—histéresis, bimodalidad, divergencia y saltos 
bruscos— se pueden observar claramente. 


Cuadro 1. Funciones elementales de catástrofe. 


Función de catástrofe Nombre 
3 A . 
x+a, x Pliegue 
4 -2 4 ¡ 
SEO nO JE Cúspide 


a atra da Cola de golondrina 


EA O TA (1 GÓ Mariposa 
a ta as ax a 
x*%y-y +0 x+0,y+a.y? Umbilica elíptica 
xy+y+0,x+a,y+a,y Umbílica hiperbólica 


xy+y+ ax+ay+a +0 Umbílica parabólica 


y2 5, ” 2, 2. 3 
xVY-Y +4,X+0,Y+AY +a,X +A. y 


2 75) 2. .2, 3 
y+y+0 x+ay+ay+a x0+ay 


ya x+ay+axy+a ya xy? 


En la figura 2 se puede ver la geometría 
de una catástrofe en pliegue, que es la más 
simple, y se empleará más tarde. En este caso, 
el conjunto bifurcación se reduce a un punto, y 
no se observa histéresis o saltos bruscos, sino 
solamente una bifurcación de las soluciones de 
la función. 


Estabilidad estructural de flujos potenciales 


La teoría de la catástrofe se ha empleado 
satisfactoriamente en varias funciones de 
diferentes tipos. Sin embargo, su aplicación 
rigurosa ha sido principalmente en el análisis de 
funciones tipo gradiente. 

Consideremos un sistema de ecuaciones f, 


que pueda derivarse de una potencial: 
0 
fi=5 vea) (7) 


donde x, representa variables de estado y a, 
son parámetros de control. La función v(x,2) 
será la función potencial. El equilibrio de ésta, 
conocido como un sistema de gradientes, 
aparece cuando: 


f=(x,0,)=0 (8) 


que corresponde a la condición de ocurrencia 
de puntos críticos, como se estableció anterior- 
mente. La teoría de la catástrofe puede aplicarse 
entonces al estudio de cambios de estado del 
sistema (los cambios en x,), en equilibrio, a 
medida que los parámetros de control a, son 
modificados. 

La función f, puede interpretarse de varias 
maneras, según el sistema real que se modele. 
Es un hecho bien conocido en mecánica de 
fluidos que si un flujo es irrotacional habrá una 
función ( tal que: 


v=grad Q (9) 


o, empleando la notación para la ecuación (7): 


(10) 
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Saltos 
bruscos 


a1 


Plano de 
control 


22 


Superficie de equilibrio 


Divergencia 


(a) 


Inaccesibilidad 


Histéresis 


(b) 


Figura 1. Geometría de catástrofe en cúspide, (a) superficie de equilibrio, 


(b) corte transversal en la superficie de equilibrio. 


donde v, indica el componente de velocidad en 
la coordenada x, mientras que el parámetro de 
control a, tendrá significado de acuerdo con la 
función potencial de que se trate. 

Ya que los componentes de velocidad, en 
términos de la función potencial, están dados 
por la ecuación (10), entonces, la condición de 
ocurrencia de puntos críticos o condición de 
equilibrio definida por la ecuación: 


v:=30(x,0,)=0 (11) 


1 


es igual a la condición de velocidad nula; esto es, 
son puntos de estancamiento. 

Se puede deducir, como resultado de la 
aplicación de la teoría de la catástrofe a flujos 
potenciales, la importancia de los puntos de 
estancamiento, porque éstos determinan la 
estabilidad estructural de los flujos modelados. 

Dependiendo del tipo de puntos de 
estancamiento de un flujo potencial, entonces 
se pueden determinar puntos de estancamiento 
degenerados o no degenerados, y el propio 
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Superficie de equilibrio 


Conjunto de singularidades 


Plano de control 


0 
1 Conjunto de singularidades 


Figura 2. Superficie de equilibrio, conjunto de 
singularidades y conjunto de bifurcaciones de la catástrofe 


en pliegue. 


flujo será estructuralmente estable o inestable, 
respectivamente. 

En el caso de flujos inestables, su función 
potencial debe poder ser reducida a una función 
de catástrofe. 


Modelo de catástrofe para un flujo 
alrededor de un cilindro circulante 


Considérese en primer término la estabilidad 
estructural de un flujo potencial clásico: el flujo 
alrededor de un cilindro con circulación, cuyo 
potencial en coordenadas cartesianas está dado 
por: 


Ua?x 
+ y 21 


0 =Ux+ (12) 


La condición para la ocurrencia de puntos 
críticos es el cumplimiento simultáneo de las 
siguientes ecuaciones: 


T T 
En 24. A Ma | 
0 >| a +1) 2 tx 


La =0 (132) 
2 
(y 
La ecuación (13.2) se cumple si: 
ax=0 (14.1) 
O si: 
0 r 2 
Y += = 14.2 
b) y/u Parar) 0 (14.2) 
que debe satisfacer simultáneamente la 


ecuación (13.1). Se puede demostrar que (14.2) 
no cumple con esta condición. 

Al sustituir (14.1) en (13.1), y después de 
realizar algunas operaciones algebraicas, se 
obtiene: 


2, P 2_ 
Uy ra =0 (15) 


cuyas soluciones quedan determinadas 


mediante la ecuación: 


1 


Dr Ey e 
a + -4 
1=2u|2n E e 


(16) 


Se obtiene así la ubicación de los puntos 
críticos, en las coordenadas dadas por: 


1 
1 T T E 22 s 

=0;y= 4u 17 

O (ES .| dd 


Las segundas derivadas parciales de la 
función potencial (ecuación (12)) son: 


IP: 
Po O SE 


7 —L, (18.1) 
(1+y?) 
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E y+ xy (sun -)-2uéx(1?- y?) 
_T 27 


% 
ay” (+?) (18.2) 
o _ 0% 
dxdy  dyox 
(24 e ylerty-y) zz (ay) 
- (+ y) (18.3) 


Es evidente, de la ecuación (17), que los 
puntos críticos ocurren en el eje y, por tanto, 
al introducir la condición x=0, las ecuaciones 
presentadas anteriormente se simplifican a: 


2 
ca E (19.1) 
dx? 
2 
00 (19.2) 
dy 
T 
| 2Ua*+— y) 
2 2 
0% _ 21 (19.3) 
dydx  dxdy al 


La matriz de estabilidad estará dada entonces 
por: 


H6= y (20) 


y su determinante igualado a cero resulta: 


2 
(2ue +35) 
21 =0 


det Hb = ra 


Q1) 


Ocurrirán puntos críticos degenerados 
cuando las ecuaciones (17) y (21) se satisfagan 
simultáneamente. 


Si se define un parámetro adimensional: 


T 


c= Ta (22) 


que equivale a decir que la circulación alrededor 
del cilindro es P=cUa, la ecuación (17) se puede 
reescribir como: 


1 
al (23) 
i 4 2|L2r 


La figura 3 muestra la gráfica para esta 
ecuación y la solución para la ecuación (21), 
para un cilindro con un radio unitario, en 
términos del parámetro adimensional c. Puede 
verse que existe solamente un punto crítico 
degenerado, el cual ocurre cuando c=4k. 

Para esta ecuación no existe una solución 
real para valores de c<4rr. Existen dos soluciones 
posibles para los valores de c>4r, de manera 
que el punto c=4rt es un punto de bifurcación. 

En la figura 4 se muestra la configuración 
del flujo para los diferentes valores de c. 


Li ds ar a a Puperficie del 
cilindro 


——- Puntos críticos 
Po... Det H/0 


Figura 3. Puntos críticos alrededor de un cilindro 


con circulación. 
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c=3571 


c=4rT 


Figura 4. Flujo alrededor de un cilindro para diferentes 


valores de C. 


Dado que existe crítico 


degenerado, el comportamiento del potencial 


un punto 


debería describirse mediante una función 
de catástrofe. Considérese una vez más la 
condición de ocurrencia de puntos críticos en 
términos del parámetro c, dada por la ecuación 
(23). Si ambos términos se elevan al cuadrado 
y se colocan a la izquierda de la igualdad, se 
obtiene: 


1712 
2 ca a C S 2 
42H H|—]-4 pa 
Y dr ¿l) E 
y si se denota: 
2 
ca a C E 2 
= + 4 = 2 
Ñ 41 2 l) e 


entonces, la ecuación (24) se identifica fácilmente 
como la ecuación de la superficie de equilibrio 
de una catástrofe en pliegue. 
Vo=y"+m=0 (26) 
Queda demostrado que el comportamiento 
del flujo alrededor de un cilindro con circulación 
se modela mediante una función de catástrofe, 
tipo pliegue, que es la más simple de todas, 
como se puede apreciar en la figura 2. 


Modelo de catástrofe en cúspide para 
flujo alrededor de cilindros y esferas 


Ya que el flujo alrededor de cilindros circulares 
es una función de catástrofe, la secuencia lógica 
es probar que lo mismo sea cierto para el flujo 
alrededor de esferas. En este caso, sin embargo, 
es conveniente usar un enfoque más cualitativo, 
usando alguna geometría de catástrofe y 
los abundantes datos experimentales que se 
encuentran en la literatura. Antes de proceder 
de esta manera es necesario mencionar 
algunos fundamentos de las convenciones de 
catástrofe. 

Dependiendo del sistema físico a ser 
analizado, es necesario establecer convenciones 
para la forma en que se lleva a cabo la transición 
entre estados de equilibrio. Estas convenciones 
se conocen como convenciones de catástrofe, y son 
principalmente dos: 


e Convención deretraso: el sistema permanecerá 
enel estado de equilibrio en que se encuentra 
(el mínimo en que se encuentra) hasta que 
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ese estado (ese mínimo) desaparezca. La 
resultante geometría de la catástrofe exhibe 
una zona de histéresis, como se muestra en 
la figura 1. 

e Convención de Maxwell: el sistema permane- 
cerá en el estado de equilibrio que minimice 
globalmente la función (en el mínimo 
global). En la resultante geometría de la 
catástrofe, la zona de histéresis desaparece, 
permaneciendo sólo la ocurrencia de saltos 
bruscos. Así lo muestra la figura 5. 


Para elegiruna convención u otra serequiere 
contar con información del fenómeno real que 
se analiza. Debido a esto, las convenciones de 
catástrofe son en realidad una extensión de 
la teoría, con miras a su aplicación a sistemas 
reales. El empleo de una u otra convención 
depende del comportamiento del sistema real. 
Si en éste no se observa histéresis, no podrá 
usarse la convención de retraso; sin embargo, 
es factible modelarlo mediante la teoría de la 
catástrofe si existen saltos bruscos, usando la 
Convención de Maxwell. 

En el flujo alrededor de objetos sumergidos, 
la variación del coeficiente de arrastre exhibe 
saltos bruscos, en relación con el número de 
Reynolds. El coeficiente de arrastre total C,, 
depende de la geometría del objeto y si no hay 
efectos de superficie libre o de compresibilidad, 
solamente del número de Reynolds: Re=VL/v, 
donde L es una longitud característica, tal como 
el diámetro de cilindros o esferas. 


EAT 


Figura 5. Geometría de la catástrofe en cúspide, 


Convención de Maxwell. 


Figura 6. Relación C,—Re para un cilindro. 


Para cada uno de los flujos, la relación 
C, Re se obtiene en forma experimental. La 
figura 6 reproduce una gráfica que determina 
esta relación para el caso de un cilindro 
circular. Se puede observar que para un Re de 
aproximadamente 4x10* existe una variación 
brusca en el coeficiente de arrastre. 

Según investigaciones anteriores realizadas 
por el autor (Martínez-Austria, 1988), en 
la aplicación de la teoría de la catástrofe, el 
segundo parámetro de control suele ser un 
parámetro geométrico, así que se procedió a 
estudiar el efecto de la rugosidad. 

Se ha demostrado (í.e. Daily y Harleman, 
1975) que la reducción brusca del arrastre en 
cilindros se verifica debido a la transición de la 
capa límite de laminar a turbulenta, de manera 
que la ocurrencia de este fenómeno depende 
de la rugosidad del objeto, lo que confirma la 
pretensión de un segundo parámetro de control 
de tipo geométrico. La gráfica mencionada 
anteriormente es válida para cilindros lisos, 
pero se han realizado experimentos para 
establecer el efecto de esta variable en cilindros 
con rugosidad. La figura 7 reproduce los 
resultados obtenidos por Fage y Warsap (1930). 
Existen otras investigaciones que han empleado 
cilindros con rugosidad distribuida (Achenbach 
y Heinecke, 1981) con resultados similares. 

Para el caso de flujo alrededor de esferas, en 
la figura 8 se muestran los resultados obtenidos 
por Achenbach (1972, 1973) para esferas lisas 
y rugosas. La reducción brusca observada en 
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o Cilindro liso s K-=7x103 
o É-=2x107 E 
D . Dm 9x107 
KE 3 
0 o LL 2x10* 


Figura 7. Relación C, Re para un cilindro con rugosidad 
(Fage y Warsap, 1930). 
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cs Liso x-5=25x10 
Y E 150x108 

o -S= 250x107 

a 5= 500x107 

e =1250x10* 


Figura 8. Relación C,, Re para esferas con rugosidad 
(Achenbach y Heinecke, 1981). 


el coeficiente de arrastre se puede interpretar 
como un indicador de catástrofe. 

La dificultad para hacer mediciones 
precisas en la zona de transición es una de las 
razones por las que la existencia de histéresis, 
por ejemplo, no se ha reportado. Por esta razón 
puede elaborarse la hipótesis de inexistencia 


de histéresis y bimodalidad, con lo que 


permanece sólo un indicador de catástrofe 
(saltos bruscos). La relación funcional entre 
variables de flujo podría escribirse de la 
siguiente manera: 

F(Cp,Re,x/D)=0 (27) 
esto es, el coeficiente de arrastre se puede 
considerar como una variable de estado, y el 
número de Reynolds y la rugosidad relativa 
como parámetros de control. 

A la luz de estas consideraciones, y de 
acuerdo con los resultados de los experimentos 
mencionados anteriormente, 
Convención de Maxwell, y los resultados de 


se adopta la 


Fage y Warsap se pueden reinterpretar para 
formar una superficie de catástrofe, como se 
muestra en la figura 9, válida para cilindros 
circulares. De manera similar, con los resultados 
experimentales de Achenbach se obtiene la 
superficie de catástrofe para esferas mostrada 
en la figura 10. 

Desde un punto de vista práctico, 
estas gráficas serían lo suficientemente 
precisas; sin embargo, debería ser posible 
establecer también la existencia de histéresis 
y bimodalidad. Shewe (1983) demostró la 
ocurrencia de estas características. Dicho autor 
empleó un túnel de viento presurizado, en el 
Góttingen Institut fúr Aeroelastik, Alemania, 
cuyas características le permitieron hacer 
observaciones de alta precisión. En primer 
lugar, debido a las dimensiones del túnel 
(60x60 cm en la sección de prueba), le fue posible 
trabajar con relaciones de bloqueo al interior 
del túnel de sólo 10%. En segundo término, 
debido a que contaba con la posibilidad de 
presurizar el túnel, le fue factible modificar el 
número de Reynolds sin afectar la velocidad 
del flujo, lo cual aumenta considerablemente 
las posibilidades de controlar la transición 
de la capa límite, al introducir un mínimo de 
perturbaciones. 

Para medir presiones, incluyendo fluc- 
tuaciones inducidas en los cilindros, Shewe 
empleó una piezobalanza construida con celdas 


piezoeléctricas, conectada a amplificadores, 
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Figura 9. Relación C,—Re representada como una superficie de catástrofe (Convención de Maxwell), cilindros. 


Figura 10. Relación C,—Re-ks/d para esferas dibujadas como una superficie de catástrofe, Convención de Maxwell. 
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para evaluar la señal en un analizador y en una 
computadora digital. 

La figura 11 reproduce algunos de los 
resultados obtenidos por Shewe. Se pueden 
observar dos zonas de histéresis y de 
bimodalidad, en el rango 10%<Re<10% que 
corresponde a la región de transición de la 
capa laminar límite a la capa turbulenta. Este 
comportamiento se observa en los coeficientes 
de arrastre y de apoyo, y en el número de 
Strouhal. 

La existencia de dos rangos de histéresis, 
saltos bruscos y bimodalidad, se puede explicar 
fácilmente si se considera que en el primer salto 
se ha producido una separación turbulenta 
sólo en un lado del cilindro; y cuando ocurre 
la separación en el segundo, se obtendrá un 
segundo salto brusco, con dos puntos de 
separación turbulenta. 

Shewe identifica como causa del compor- 
tamiento del fenómeno la inestabilidad 


Sp 0.4 j Ml 


¿tá—Á + 


0.5 


105 106 Re 107 


Figura 11. Variación del coeficiente de arrastre, número de 
Strohual y media cuadrática del coeficiente de apoyo para 
un cilindro, en la zona de transición (Shewe, 1983). 


hidrodinámica y lo interpreta como transiciones 
de fase entre los dos estados de orden, 
proponiendo diagramas de bifurcación. 

Debe señalarse que estos rangos de histéresis 
son muy pequeños, incluso en condiciones 
experimentales bien controladas; de manera 
que en términos prácticos es conveniente 
emplear la Convención de Maxwell, como se 
mostró anteriormente. 


Conclusiones 


Mediante la aplicación de la teoría de la 
catástrofe a los flujos irrotacionales, se ha 
demostrado analíticamente que el flujo alre- 
dedor de cilindros circulares con circulación es 
una catástrofe en pliegue. 

Este resultado se ha ampliado, proponiendo 
que el flujo rotacional alrededor de cilindros y 
esferas también se puede modelar mediante la 
teoría de la catástrofe. Para ello, la rugosidad 
relativa se eligió como un segundo parámetro de 
control y, empleando resultados experimentales 
disponibles en la literatura, se demostró que 
ambos flujos pueden representarse como una 
geometría de catástrofe de tipo cúspide. Debido 
al pequeño rango de histéresis en estos flujos, se 
propuso como más conveniente la Convención 
de Maxwell para estos casos. 

Con ambos resultados se puede formular 
la hipótesis de que los flujos alrededor de 
objetos sumergidos constituyen procesos de 
catástrofe y que este tipo de comportamiento 
debiera extenderse a cualquier tipo de 
geometría. 


Nomenclatura 

a Radio de cilindro. 

a, Parámetros de control para una 
familia de funciones; j=1, 2, 3,..., k. 

b,c,d Constantes. 

Es Coeficiente de arrastre total. 

Co Coeficiente de presión base. 

d Ancho de estela. 

det Determinante de la matriz. 

grad Gradiente de una función. 
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Función de n variables. 

Familia de funciones con n 
variables y k parámetros. 

Matriz de estabilidad o hessiana 
de la función f y 6. 

Longitud característica. 
Exponentes y constantes. 
Rotacional de una función. 
Número de Reynolds. 

Número de Strouhal. 

Velocidad de flujo uniforme. 
Velocidad media de un flujo. 
Componentes del vector 

v, en relación con x, y, z, 
respectivamente. 

Variables de una familia de 
funciones; ¡=1, 2, 3,... 1. 
Coordenadas en un sistema de ejes 
tridimensional. 

Vector de velocidad con 
componentes u, 0, Ww ÍN X, Y, Z, 
respectivamente. 

Función potencial, de manera tal 
que v=grad d. 

Circulación. 

Rugosidad relativa. 

Coeficiente de una forma 
cuadrática, autovalores 
(eigenvalores) de la matriz de 
estabilidad. 

Viscosidad cinemática del fluido. 
Punto crítico de una función. 
Operador matemático. 

Denota “evaluado en el punto p”. 
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Abstract 


MARTÍNEZ-AUSTRIA, PF. Catastrophe models for the flow around immersed bodies. Water 
Technology and Sciences, formerly Hydraulic Engineering in Mexico (in Spanish). Vol. 
L no. 1, January-March, 2010, pp. 7-19. 


In fluid mechanics there are several phenomena that exhibit a hysteretic behavior and sudden 
changes not only of a quantitative, but also a qualitative nature. Such phenomena can be 
represented and explained by the catastrophe theory, which explains this kind of behavior 
observed in some physical systems. This article offers an analytical demonstration that the 
flow around immersed cylinders is a fold catastrophe. Extending this result, now with a semi- 
empirical approach, it can be shown that flow around cylinders and spheres with roughness 
can be represented using a cusp catastrophe geometry. Both results support the hypothesis that 
the flow around objects immersed in liquids can be modeled through catastrophe geometries 
and hence explain the occurrence of bifurcation, bimodality, sudden jumps and hysteresis. 


Keywords: catastrophe theory applications, catastrophe theory, flow around immersed bodies, 
hysteresis in fluid mechanics, fluid mechanics, cusp catastrophe. 
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